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Abstract 

The paper deals with the study of superfluidity by a Ginzburg-Landau 
model that investigates the material by a second order phase transition, 
in which any particle has simultaneouly a normal and superfluid motion. 
This pattern is able to describe the classical effects of superfluidity as the 
phase diagram, the vortices, the second sound and the thermomechanical 
effect. Finally, the vorticities and turbulence are described by an extension 
of the model in which the material time derivative is used. 
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1 Introduction 

The behavior of superfluids is very different from the phenomenology of a perfect 
fluid, as well as superconductivity [14] is a phenomenon very different from per- 
fect conductivity (see Bardeen [4], Chandrasekhar [6], Landau [17], [16], London 
[19], Mendelsshon [20], Fabrizio,Gentili and Lazzari [10]). Meanwhile, there is 
an evident similarity between the behavior of superfluids and superconductors. 
By the way, London claims in [19] § 22 "... in either case, superconductivity and 
superfluid helium, the basic equations can be written in the same form" , and 
Mendelsshon asserts "However, frictionless flow and persistent currents are not 
the only features in which He H resembles a superconductor. As at the transi- 
tion of a metal to the superconductive state, there occurs in helium also a rapid 
drop in the entropy. Moreover, superfluid flow, just as a persistent supercurrent, 
is distinguished by zero entropy" . In view of the apparent analogies between 
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superfluidity and superconductivity, in this paper we present a phase transition 
model, which use the two fluids Ginzburg-Landau view point, in which the equa- 
tion governing the motion of superfluid component Vs assumes an analogous 
form as the equation for superconducting electrons. It turns out in our analysis 
that V X Vg ^ 0. Accordingly, the equations allow to predict existence of vor- 
tices, which have been widely studied by several authors (see: [21] cap. 6,[1], 
[2], [3], [5]). Moreover in this paper, we characterize superfluidity as a second 
order phase transition^. Thus, our model for superfluidity is able to explain the 
phase diagram represented in Fig.l, at least in the region close to the A— line, 
which is the most significant area. At the same time, it allows to prove the 
existence of a critical velocity of superflow, beyond which the frictionless flow 
breaks down. 
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Fig.l. Helium II. Phase diagram. 

A remarkable feature of our model is the constraint on the rotational motion by 
means of the equation V x — v„ which relates Vs with the normal component 
of the velocity v„ . So that the velocity is defined by 

V = v„ f^r,. (1) 

The formula (1) does not allow to define the superfiuids by means of two compo- 
nents, "the treatment of the liquid as a mixture of normal and superfiuid parts 
is simply a form of words ....Like any description of quantum effects in classical 
terms, it is not entirely adequate. It does not at all mean that the liquid can 

^ There is a marked difference between first-order and second-order phase transitions. In 
first-order transitions, below the critical temperature the whole body occurs in a single phase. 
Instead, in second-order transitions, below the critical temperature the body consists of the 
simultaneous occurrence of both phases. That is why we can view the material, below the 
critical temperature, as consisting of both the normal fluid and the superfluid (two fluid 
model). 
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actually be separated into parts" [15]. Actually as in (1), any superfluid particle 
is capable of two types of motion simultaneously, corresponding to the energy 
spectrum of phonons and rotons. "One of these motions is normal, i.e. has the 
same properties as that of an ordinary viscous liquid; the other is superfluid" 
[15]. 

In Section 2 we will examine thermodynamical consistence of the model. 

Then, in view of the similarity between superfluidity and superconductivity, we 
will prove existence of vortices, whenever a superfluid, contained in a cylindrical 
bucket, rotates rapidly around its axis. In particular, in accordance with exper- 
imental results, we will show that, if the angular velocity is low, the superfluid 
stays at rest. Instead, when the velocity exceeds a threshold value, there occur 
vortices, similar to the vortices observed in superconductivity. Moreover, the 
proposed model allows to explain thermomechanical effect, in which we observe 
a motion of liquid produced by heat, but in the same direction to heat flux. 
Finally, we prove the existence of the phenomenon of second sound. 

In section 3, we investigate the problem which arising when the flow is such 
that we need to consider the material time derivative of the velocities. This 
study introduces new non-linear terms that increase the complexity of the model 
and require a suitable change of the previous pattern. 

In the last part of the paper we study the connection between this model and 
the turbulence. So that, as observed by Mendelsshon [20] "It is now clear that 
the dependence of the heat conduction on the heat current originally observed in 

Cambridge in 1937 is evidence of turbulence" . A remarkable feature of our 

model is the constraint on the rotational motion V x = v„, that relates 
with the normal component of the velocity v„ and which provides the behavior 
of a disordered set of thin vortices. 



2 A first model for superfluidity 

In this section we shall provide a first approximate model for the study of 
superfluidity in a domain O c as a second order phase transition, by means 
of the Ginzburg-Landau equation 

§ = -VV-/(/^-l+« + Ap + vf) (2) 
at K 

where / denotes the order parameter (or phase field), while p is the pressure, u is 
the absolute temperature and k, A are positive constants. As studied in [8], the 
Ginzburg-Landau equation has to be considered as a new field equation, that 
we obtain by the balance law on the structure order (See Appendix). Indeed, 
following Landau, we suppose that the transition of the Hell from a normal to 
a superfluid state induces a change in the internal structure order. The phase of 
this transition is represented by the scalar parameter / e [—1,1] that is linked 
with the density of the superfluid particles by the formula 

= f. 
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Hence the phase field = denotes the normal state, while G (0, 1] describes 
a superfluid state. In this framework, in equation (2) if 



R = u + Xp + v;-l>0, 

then the fluid is in the normal state, while when R <0 the pattern describes the 
superfluid phase. In this work we do not employ the classical two fluids model 
which introduces two different velocities and densities, as for a mixture, but we 

suppose the velocity of any particle composed of two excitements, the normal 
component v„ and the superfluid component v^, such that 

v = v„ + /2v«. (3) 

Because superfluidity as well as superconductivity must be studied as a second 
order phase transition, then under the transition temperature, the phase of some 
particles / can be still equal to zero. Nevertheless, according as the particles 
are in the normal phase (/^ = 0) or in the superfluid one (/^ > ), we cannot 
interpret this different behavior as two different fluids. 

Henceforth in this paper we choose A = 0, and the density p a positive 
constant. In addition, in this first approximate model wc suppose the motion 
such that the acceleration a„ = Then, for the velocity v„ we propose the 
modified Navier- Stokes equation 

= -Vp - mV X V X v„ - mV X fv, + g (4) 

where g denotes the external force and fj, is the viscosity coefficient. Since we 
have supposed the superfluid as an incompressible material, then the continuity 
equation provides 

V • v„ = (5) 
Furthermore, we suggest the following equations for the component 

= - X V X - jifws + Vu + Mh (6) 

V • f\s = -rf^ (7) 

where ^ is a suitable scalar function and r is a positive constant. While the 
vector h is such that 

/xV X h = g - VjJ , V • h =0, h X n\g^ = (8) 

Now, using (8), it is easy to prove that equation (4) can be obtained from 
equation (6) if 

V X v, = v„ (9) 

Therefore, under the latter condition it will be sufficient to work only with 
system (6)-(9), which contains equation (4) as a consequence. 
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3 Thermodynamics 



In order to obtain the heat equation, let us consider the first law of thermody- 
namics in the form 

E = V^ +V^'' + h, (10) 

where E is the total energy and h is the rate at which the heat is absorbed by 
the material. The internal powers V-^^V^", related to the variables /, v^, are 
given^ by 

r^' = f! + j-J{^f)\ + \[ii-m, + fMu + vl) (11) 

From first law (10) we have 
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Fig. 2. Specific heat of Helium. 



If T = iv^ + i/^ denotes the kinetic energy, and 

e = J c{u)du + \{l f f + , (14) 

is the internal energy, then the total energy E ^ e + T is given by 

E = + If + / c(u)du + 1(1 - ff + ^(V/)^ , (15) 

■^See Appendix. 
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In our model the Fourier theory of heat conduction wiU not be modified. Ac- 
cordingly, the constitutive equation relating the heat flux q to the gradient of 
the temperature, assumes the classical form 

q = -k{u)Vu, (16) 

where the conductivity k depends on the absolute temperature. In this frame- 
work, since the conductivity is very small when the absolute temperature is 
close to zero, we suppose that 

k{u) = kou , fco > . (17) 

The heat balance is expressed by the equation 

h= -y ■ q + r, (18) 

from which we deduce the entropy equation (see [7]) 

--^ q-Vw = -V--5 + -. (19) 

Therefore, 

= - - 4 q • Vu (20) 

can be identified with the rate (per unit volume) at which entropy is absorbed, 
while q/u and r/u represent the entropy flux and the entropy supply. Then 
equation (19) describes the entropy balance [7], [11]. Because of the non-local 
constitutive model given by system (2-7), we can state the following (see [7], 
[9]) 

Second law of thermodynamics (non-local form). There exists a function 
r]{u,f), called entropy function, such that 

5^>ii' + V- k= -- 4q-Vw+V-k. (21) 

at u 

where k is a suitable flux vector such that k • nj^j^ = 0. 

It is easy to prove that system (2-7) with constitutive equations (13), (16) 
agrees with the inequality (21), if we regard as entropy function the functional 

V{u,f)= j^du-G{f). (22) 

where G{f) = ^/^and as flux vector k =f'^Vs- Finally, in the sequel we shall 
consider an approximation of (19), namely 

- G{f) = V • (/2 V,) + fcoV^w + - , (23) 

u u 

where we have neglected the dissipative terms f^, (vgt+Vc/)— Vw)^, {WuY ^rf^cj)^ . 
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Now, let us evaluate the free energy defined by 

tjj = e — UT] . 

From expressions (14) and (22) of the internal energy and the entropy, it follows 
that 

i>=^{Vfr + \{l-fr+uG{f) + {l-u) j ^du. (24) 

We will prove that system (2), (3), (4), (23) allows to describe the main aspects 
of superfluidity. First, we examine the phase diagram represented in fig.l. If 
we focus our attention on phenomena in which the pressure is below 25 atm, 
the transition occurs along the A— line, which we can assume close to a vertical 
line. This phenomenon is explained by equation (2). Indeed, by determining the 
transition states {ut, v^t), we observe that the line which separates the normal 
and the superfluid phases is given by the curve 

u + Xp + vl = l (25) 

If we suppose A <C 1 and consider the equilibrium states, i.e. Vg = 0, then from 
(25) we get the equation u = I, which is a good approximation of the A— line 
(see fig.l). Instead, when ^ 0, the model is naturally able to account for 
the existence of a critical velocity, namely a threshold value above which the 
frictionless flow breaks down. This follow directly from equation (25). 



The latent heat I related to the transition is given by 

I = Ut{v{ut, fi) - v{ut, /2)) 

where /i and /2 are the local minima of the free energy at the transition points 
u + Vg = 1, In this case fi — f2 = 0. Hence, I = 0, which means that the 
transition is of second order. 

In summary, our differential system assumes the form 

^^Iv^f-fif-l + u + vl) (26) 

= - X V X vg - fj,fvs + S/U + iih (27) 
V • /'v, = -r/2<^ (28) 

^|-/| = V.(/^v.) + fcoV^.+ ^ (29) 
u at at u 

Moreover, we associate to system (26-29) the boundary conditions on the domain 
O 

V/ • n|an = , (V x v^) x n|gn = w (30) 
^ 0> v„|gf^ = 0, u\Q^ = Ub (31) 
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It is worth noting that, when / = 0, in view of (3), v = v„ and equation 
(6) is inessential, since the motion is described by means of the Navier-Stokes 
equations 

dv 

— = -Vp -/iVxVxv + g, V-v = 0, (32) 

obtained from (4), (5). On the contrary, when the transition occurs, / > and 
the consequent breaking of symmetry leads to system (26)- (29). Well-posedncss 
of differential problem (26)- (29) in the case / = has to be investigated carefully 
in order to prove uniqueness of solutions. Indeed, when / = equation (27) 
reads 

= - /xV X V X - Vu + Mh, (33) 

while (28) is satisfied identically. The temperature u is univocally determined 
through (29), endowed with suitable initial and boundary conditions. However, 
since equation (33) is no longer coupled with (28), uniqueness of solutions cannot 
be proved. For this reason, it is necessary a choice of a gauge, which consists in 
decomposing the velocity Vg as 

Vg = A + Wif, or V X A = v„ 

The vector A is specified by fixing a value for its divergence and choosing suit- 
able boundary conditions. By analogy with superconductivity, we will choose 
London's gauge 

V- A = 0, A - n|an = 0. 
Therefore, from (33) we deduce 
dA 

—= - V((/> + (f) - IJ.V X \/ X A - \/u + iJ.h (34) 
Finally, from equation (34) we get equation (32). The main result which this 

model allows to prove is the presence of quantized vortices, when a cylindrical 

bucket containing Helium II rotates around its axis. Consider the stationary 
problem, in addition let us suppose the temperature gradient Vu = 0. Hence 
the differential system (26)-(28) with h = assumes the form 

-V'f~fif-l + u + ivl) = (35) 

-V(t> - IJ.V xV xvs - nfvs = (36) 
V • V, + r/2</. = (37) 

together with the boundary conditions 

V/ • n\en = , (V x Vg) x n\an = w , ■ n\an = 0. (38) 

If the angular velocity is below a critical value , depending on the geometry 
of the problem, the superfluid component stays at rest, with the exception of 
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a very small penetration depth close to the boundary^. This equilibrium state 
is called Landau's state. The phenomenon is analogous to the diamagnetism 
shown in superconductivity by the Meissner effect, when for a second type 
superconductor the magnetic field is lower that a given field iJci ■ When the 
angular velocity is greater that , we observe the creation of vortex lines in an 
exactly similar way as in superconductivity. Since, problem (35)-(38) is identical 
to the Ginzburg-Landau system in stationary case, studied by Abrikosov [1]. for 
superconductors of type II. Therefore, when uj > cOci with the same analysis, we 
may prove the existence of vortex lines in a superfluid. Moreover, as observed 
in [21], when the rotation velocity exceeds a second threshold value ujc2, owing 
to increase of vortices and to the overlapping of their cores, superfluidity is 
destroyed. Even this effect is analogous with the corresponding phenomena in 
superconductivity and it can be explained by (34). Indeed, for sufficiently large 
values of Vg we have 

u+ v2 > 1 , 

so that the fluid is in the normal phase. In every practical situations, angular 
velocities u satisfy the condition 

CO < Wc2 ■ 

Moreover, the flrst critical velocity is such that 

Wcl < 1 . 

Another important phenomenon occurring in superfluidity, it is the thermome- 
chanical effect according to which the particles of the superfluid flow in the 
same direction to heat flux. This phenomenon, which contradicts the classical 
behavior of thermo-fluid mechanics, can be explained by means of equation (6). 
Indeed in (6) the sign of Vu is the same as the acceleration Vgt of the superfluid 



^ Consider a fluid in a long straight wire of circular cross-section carrying a rotation around 
the axis. We suppose cylindrical symmetry and indroduce the coordinates z, r, (k, r, j). Then 

we have 

Vs{r) = ve{r)i 
4> = Hr) 



v„(r) = Vxv.(r) = ^^^k 
or 



Therefore, from (36) we obtain 

V</)(r) = 

-72- = — :r'^^f 

ar-' r dr 

The solutions related with a given fleld vq at the surface of the cylinder of radius R is 

h{r) 



Vg = Vo 



h(R) 



where /i is the modified Bessel functions of imaginary argument. Because of the behavior of 
Ii (r) for r < R, the velocity vg (r) will be confined to a thin layer next to the boundary. 
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component. This means that the gradient of the temperature causes an increase 
of the velocity in the same direction. Moreover, to understand better the 
phenomenon, we have to observe that in the thermo-mechanical effect the tube 
may be straight and very narrow. As a consequence, the motion of the super- 
fluid phase (/ ^ 0) through the tube is required to satisfy Vx = v„ = 
and hence the velocity is given by v = /^Vg. In other words, since the normal 
component is absent, the motion is necessarily that of the superfluid component. 
This in turn allows for the flow through very narrow tubes. It follows from (6) 
that, in stationary conditions, the motion is governed by the equation 

-V<^ - f^s + Vu + h = 0, (39) 

which, consistent by (39), shows that the term Vm favors the particle dis- 
placement toward regions at a higher temperature. Moreover, because of the 
narrowness of the tube, only superfluid particles (/ — 0) are allowed to flow in. 
Instead the normal component undergoes to a viscous resistance which forbids 
the flow to cross the tube. 

Finally, we study the second sound effect. For this problem, let us suppose 
/ = const, and the supplies h = 0, r — 0. Then, from (27)-(29) it follows that 

V • v«t = -V20 - /iV • fvs + V^w (40) 

autt = V ■ (/' v,t) + koV^ut (41) 

where a = will be supposed a positive constant. Moreover, in the superfluid 
framework we can suppose the thermal conductivity ko ^ 1, therefore we get 
from (37) and (41) 

In other words, the temperature satisfies a wave equation that is able to provide 
a new kind of wave. 



4 General case 

In the previous sections the material time derivative ^ = ^ + v • VF was 
always completely ignored, because we have supposed slow motions. In this sec- 
tion we shall study the general case and then we use the material time derivative 
^ which we denote only with F. In addition we suppose the fluid incompress- 
ible with a given constant density p = 1. Then, on the domain O, the differential 
system assumes the new form 

/= V/ -/(/'-! + « + v,^) (42) 

v„ = -Vp - AtV X V X v„ - /iV X fvs + g (43) 
V • v„ = 
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= -V0 - /xV X v„ - nfvs + Vu + /ih (44) 

V • fv, = -Tfcf> (45) 

Becaiisc of the use of the material time derivative, even if we suppose /iV x 
h = g — Vp, it is not possible to prove that from equations (43) and (44) we 
obtain equation (4). Otherwise, if we consider the curl of the equation (44) and 
compare this new equation with (43), we have the restriction 

V X Vs = v„ (46) 

Moreover, the internal powers , V^^, assume the analogous form 

V' = + ^ [(V/)^]' + J [(l-ffl^ +/./(" + V?) (47) 
V^' = /iVx v„-v, + ^-i(v, + V0-Vu)2+/2v, -v, +mV-/2v, +r 
In particular, from (43-44) we have 

= b'l] ■ + m(V X v„ + /2v,) • V X v„ - (g - Vp) • v„)da; 



= Inih [v^] ■ - V X v„ • v«)rfa 



(49) 



Hence, if (9) holds the identity (49) is satisfied, but the vice versa is not true. 
However, from (49) and (31) we obtain the restriction 

_£(^[vn]-Vxv,-v„)dx = (50) 

which is satisfied by means of (46). 
Then, from the first law (10), we get 

^ = ^ + /' + ^[(W)']' + ^[(l-ffl'+//(« + v2)+Vxv„.v,(-^l) 

where the internal energy is given by 

e = j c{u)du + 1(1 - ff + i-(V/)2 , (52) 

which implies 

h = c{u)u -P- (v, + V<P- Vuf - uff - uV ■ fv, - Tpcf , (53) 

where c{u) = eu{u) is the specific heat. 

With similar arguments considered in the previous section we obtain the 
same function for the entropy 



V{uj)= I ^du-G{f). 



(54) 
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and under the same approximations, the following heat equation 

_ G(/) = V • (/2 v,0 + fcoV^u + ^ , (55) 

Prom expressions (15) and (54) of the internal energy and the entropy, it follows 
that 

i>=^{Vff + \{l-ff+uGif) + {l-u) I ^du. (56) 

We will prove that system (42)-(45), (55) allows to describe the main aspects of 
superfluidity, but in more complex form, because of the presence of total time 
derivative and of equation (50). 



5 Turbulence 

It is well known that turbulence is a phenomenon which has been observed 
quite frequently in superfluidity. " Hydrodynamic flow in both classical and 
quantum fluids can be either laminar or turbulent. To describe the latter, 
vortices in turbulent flow are modelled with stable vortex fllaments. While 
this is an idealization in classical fluids, vortices are real topologically stable 
quantized objects in superfluids. Thus superfluid turbulence is thought to hold 
the key to new understanding on turbulence in general." [5], (see also [12]). 

This efl:ect is connected with the deflnition (3) of velocity and with the rela- 
tion (9) or (50) between normal and superfluid components. In particular, if we 
consider the model described by the system (2), (3), (4), (23), then we observe 
that, under the critical temperature T^, the path of the particles assumes an he- 
licoidal motion, whose circulation, defined around the vortex core, is quantized. 
If the temperature goes below a second critical temperature Tc2, the number 
of the superfluid particles increases. In such a case the approximation (2), (3), 
(4), (23) is not correct, because we need to consider the time material deriva- 
tive. Hence the motion assume a chaotic form, typical of the turbulence flow. 
Otherwise, it is possible to approximate the system (42-45), (55) considering 
the material time derivative for v„ and Vg, then we obtain the system 



^ = -VV-/(/^-l + « + v?) (57) 
at K 

Vn = -Vp - /xV X V X v„ - mV X fv, + g (58) 

V X V, = v„ , V • v„ = (59) 

= -V(/) - /LtV X v„ - m/^Vs + \/u + jih. (60) 

V • f^s = -rfcf (61) 

(ill • v 

co^-G{f)=W■{f^rs) + koV\+- (62) 



where in this framework we have supposed = cq positive constant. This 
system is completed by the boundary conditions (30-31). 
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Now, let us consider the phenomenon of turbulence in viscous fluids. As for 
a superfluid we suppose the velocity given by (3). While, we suggest for the 
differential problem a first order phase transition model by the following system 

^ = l^'^ ~ + (63) 

v„ = -Vp - AiV X V X v„ - /iV X G(/)vs + g (64) 

V X v« = v„ , V ■ v„ = (65) 

V, = -V(l> - mV X v„ - iJ,G{f)vs + Mh (66) 

V • G(/)v, = -uGif)<p (67) 

where g and h are related by (8), Nr is a constant related with the Reynolds 
number and ly a positive constant. Moreover, because of the first order phase 
transition described by the system (63-67), the functions F{f) and G{f) are 
now given by the following fourth order polynomials 

fi f?, fi nfS f2 

Fit) = — - — , G(/) ^- ^—^L_ (68) 

Finally, in this framework the internal powers connected with the variables 
/, and Vs are given by 

= f' + ^[{^ff]'+NRF{f)+G{f)vl (69) 
V-^ = G(/)v,.v,+/x-i(v. + V</))^ + ^[v^]- (70) 
= ^(Vxv„)2 + G(/)v,.Vxv„ + l[v2]- (71) 

For this isothermal problem the laws of thermodynamics are given by the 
following 

Dissipation Principle. There exists a functional '^(/, Vg), called free en- 
ergy, such that 

V' < + ^ [(V/)^] • + NnF{f) + n-\v,+Vcl>f 

Then, the free energy is function only of / and is defined by 

V'(/) = ^(V/)^ + 7VflF(/) 
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6 Appendix 



In many natural phenomena, as phase transitions, chemical reactions, biological 
(processes) models, we can observe a natural evolution toward change of ma- 
terial structure order. We suggest that this behavior is related with a natural 
process which we may represent as a balance law on the structure order and it 
will be given as a function of order parameter / . Moreover, as for the other field 
equations, to this new law we may connect an internal power, which must be 
considered in the energy balance law. In this framework the structure order is 
a new form of energy, because during the transformation we observe a variation 
of the structural energy. Actually, if we consider a phase transition we no- 
tice a transformation from a less ordered material structure to a more ordered 
one or vice versa (see [8], [13]). Moreover, below a critical temperature, the 
structure order of many materials is greater then above. We meet, analogous 
behaviors, during biological precesses and chemical reactions, always connected 
with structure order variations. 

In order to obtain a balance law on the structure order. Consider a body 
B, for any sub-body C 6, we denote with S^{S) the rate at which structure 
order is absorbed by the material per unit time, given by 

S\S) = [ pkdv (72) 
Js 

where p is the density and k the internal specific structure order. While, the 
external order structure S^{S) assumes the form 

5^(5) = / p-nds+ [ pSdv (73) 

JdS Js 

where the vector p denotes the order structure flux and 6 the structure order 
supply. 

Hence, the order structure balance is given for any S C B hy the equality 
/ pkdv = I p ■ nds + / p5dv (74) 

Js JdS Js 

In local form the integral equality (74) implies the identity 

pk = V ■p + pS (75) 
In any model of phase transitions the functions k and p are usually defined by 

k = ft + F'{f)+mG'{f) (76) 

P = -V/ (77) 

K 

where m is a suitable coefficient depending of the field. While F and G are 
functions that characterize the order and the feature of the transition. For a 
second order phase transition as superfiuidity F and G are defined by 
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while m = u+Ap+Vg, from which we obtain by (75), (76) and (77) the equation 

Pft = -V^/ - pF'iJ) - p{u + \p + ^l)G'{f) (79) 

i.e. the equation (2) with p = 1 

Because the equation (75) is joined with an energy, we have to introduce the 
power balance connected with this equation, namely 

pkft + p • V/t = V • {ftv)+pftS (80) 

Finally, we denote by 

= pkft + p-Vft, vi = V-{ftp)+pftS (81) 

the internal and the external order structure power density respectively. 
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